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Abstract—The letter is devoted to the introduction of a new tech- 
nique to derive the cloud bottom height (CBH) from satellite mea- 
surements of the cloud reflectance in the oxygen A-band. The infor- 
mation on the cloud top height needed in the retrieval of the CBH 
must be obtained from separate measurements to insure small bi- 
ases in the retrieved CBH. Such measurements can be performed 
by a space-based lidar. 


Index Terms—Cloud remote sensing, O2 A-band spectrometry, 
radiative transfer. 


I. INTRODUCTION 


HE DETERMINATION of cloud geometrical character- 

istics (e.g., the cloud geometrical thickness /, the cloud 
top height A, and the cloud bottom height H) is of a great im- 
portance for climate issues, dynamical meteorology, and air- 
craft safety issues. Global determinations of geophysical pa- 
rameters are possible only using satellite measurements. There 
are a number of techniques to determine the cloud top height 
from a satellite. A recent review of these techniques is given 
by Rozanov and Kokhanovsky [12]. The techniques range from 
active measurements (e.g., see Winker et al. [15] and Poole et 
al. [9]), using lidars and radars on satellite platforms, to passive 
techniques, based on the processing of information, contained 
in the thermal infrared radiances [4], Ring effect [2], reflected 
light polarization [5], [14], stereo photogrammetric technique 
[8], and oxygen A-band technique [12], [16], to name a few. 

Rozanov and Kokhanovsky [12] has extended the Og A-band 
technique to the determination of the cloud geometrical thick- 
ness. However, it was found that simultaneous determination of 
both the cloud top height and the cloud geometrical thickness 
introduces large biases in the retrieved cloud bottom height in 
many cases. 

To reduce these biases we propose to use the simultaneous 
measurements of the cloud top height by a lidar system (the ac- 
curacy in h is better than 20 m) and the cloud reflectance spec- 
trum by the O> A-band spectrometer. These two independent 
pieces of information can be used to determine the cloud bottom 
height H = h — l from a satellite or aircraft. The study of the 
accuracy of such a technique is the main subject of this work. 


Il. THEORY 


The cloud reflectance spectrum is almost not sensitive to the 
cloud geometrical characteristics outside of gaseous absorption 
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bands. The situation is radically changed if we consider the ra- 
diative transfer in the molecular absorption line [12], [13], [16]. 
Indeed, let us assume that we have a gas in a planetary atmos- 
phere, which absorbs almost all incident radiation in a narrow 
band. Then the depth of this band, measured by a receiver on a 
satellite will depend on the cloud altitude. Gas concentrations 
generally decrease with the distance from the ground. There- 
fore, clouds at a high altitude do not allow most of photons to 
penetrate to low atmospheric layers and be absorbed there. So 
the depth of a molecular line in the reflected light will decrease, 
if high clouds are present in the field of view of a sensor. 

The next question to address is the influence of the cloud 
geometrical thickness on the reflectance spectrum R(A) in the 
gaseous absorption band. One expects that spectra P(A) in the 
gaseous absorption band for clouds having the same top heights 
but different cloud geometrical thicknesses will differ even if the 
cloud optical thicknesses 7 coincide [6]. This is due to the fact 
that multiple light scattering will lead to larger average photon 
path lengths in clouds as compared to a cloudless atmosphere, 
thereby increasing absorption. This must lead to the decrease 
of the reflectance for geometrically thicker clouds. Radiative 
transfer calculations (see Fig. 1) confirm this fact. 

The spectrum /() as shown in Fig. 1 can be used to deter- 
mine the cloud bottom height from a satellite. The main steps 
of the inversion technique are given below. 

First of all, the TOA reflectance R is presented in the form 
of a Taylor expansion around the assumed value of the cloud 
bottom height equal to Ho 


R(H) = R(Ho) + ` a;(H — Ho)! 


i=l 


(1) 


where a; = R“(Ho)/i!. Here R (Ho) is the i-derivative of 
R at the point Hg. The next step is the linearization, which is a 
standard technique in the inversion procedures [11]. We found 
that the function R( H) is close to a linear one in a broad interval 
of the argument change [6]. Therefore, we neglect nonlinear 
terms in (1). Then it follows that 

R = R(Ho) + R'(Ho)(H — Ho) (2) 
where R’ = (dR)/(dH). We assume that R is measured at sev- 
eral wavelengths (A1, A2,..-An) in the oxygen A-band. Then 
instead of the scalar quantity FR we can introduce the vector 
Rimes With components (R(A1), R(Ag),... R(An)). The same 
applies to other scalars in (1). 

Therefore, (2) can be written in the following vector form: 


j = ax (3) 
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reflection function 
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Fig. 1. Cloud reflection function in the oxygen A-band calculated using 
SCIATRAN (lines) [10] and asymptotic analytical theory (symbols) [6] at the 
nadir observation, the solar zenith angle equal to 60°, the optical thickness equal 
to 50, the cloud top height equal to 9 km, and the cloud bottom height equal 
to 3 km (circles) and 7 km (triangles). All other parameters (e.g., atmospheric 
vertical profiles) coincide with those described in [6]. The reflectance function 
is averaged with the step 0.2 nm using a SCIAMACHY [1] response function. 


where 7 = Rmes—R(Ho),@ = R'(Ho), and £ = H— Ho. Note 
that both measurement and model errors are contained in (2). 
The solution £ of the inverse problem is obtained by minimizing 
the following cost function: 


è = |7- del]? (4) 


where || || means the norm in the Euclid space of the correspon- 
dent dimension. 

The value of ĉ, where the function ® has a minimum can be 
presented as 


(5) 


where (,) denotes a scalar product in the Euclid space, n is 
the number of wavelengths, where the reflection function is 
measured. 

The functions R(Ho) and R’(Ho) in (5) must be calculated 
using the radiative transfer theory with input parameters char- 
acteristic for a given atmospheric state. We use the approximate 
analytical theory for such a calculation. The details of this theory 
are given by Kokhanovsky and Rozanov [6]. Basically, the ap- 
proximation has an accuracy better than 5% as compared to 
line-by-line calculations in the O2 A-band [6] for typical cloudi- 
ness with t > 5. The use of the exact radiative transfer theory 
is also possible but it leads to huge calculation time and does 
not provide a better approach to the problem at hand due to all 
uncertainties involved (e.g., possible multilayered cloudiness). 
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Therefore, knowing values of the measured spectral reflec- 
tion function Rmes and also values of the calculated reflection 
function R and its derivative R’ at H = Hp and several wave- 
lengths, the value of the cloud bottom height can be found from 
(5) and equality: H = ĉ + Ho. The value of Ho can be taken 
equal to 0.5 km, which is a typical value for low level clouds 
at middle latitudes [3]. The main assumption in our derivation 
is that the dependence of R on H can be presented by a linear 
function on the interval x [6]. 

The retrieved value of H is compared with Hp. If the differ- 
ence is smaller than 100 m, the value of H is taken as a retrieved 
value and the inversion procedure is finished. Otherwise, the re- 
trieved value of H is substituted in (2) instead of Ho and itera- 
tions are performed until the convergence is reached. We need 
three to five iterations to reach a convergence, if the initial guess 
is far away from the true value of the cloud bottom height. 

Several additional parameters are needed in the retrieval 
procedure. They are the cloud optical thickness, the cloud 
liquid water profile, the cloud thermodynamic state, the cloud 
droplet radii, etc. They must be simultaneously derived or 
assumed using climatological values. We have found that the 
cloud optical thickness 7 is the most important parameter, 
which influences the retrieval. So we find the value of 7 from 
measurements outside the gaseous absorption band as described 
by Kokhanovsky and von Hoyningen-Huene [7]. 

Clearly, if the forward and inverse models use the same 
system of equations, the inverse problem solution reproduces 
accurately the input parameters for the forward model in the 
absence of the measurement noise. We have checked this using 
the forward and inverse models based on the same set of analyt- 
ical equations as described by Rozanov and Kokhanovsky [12]. 
Indeed, the cloud bottom height used as input in the retrieval 
scheme coincided with the cloud bottom height retrieved by the 
solution of the inverse problem as specified above in this case. 

The next possible step is to introduce measurement errors and 
see the influence of these errors on the retrieval of H. How- 
ever, we have chosen a different strategy in this work. Namely, 
we calculate the cloud reflectance spectrum R(A) in the Og 
A-band using the exact radiative transfer calculations with the 
latest version of the radiative transfer solver SCIATRAN [10] 
and use this exact spectrum in the analytical retrieval procedure 
described above. Because possible measurement errors are well 
below the accuracy of analytical equations (see Fig. 1), such 
an approach can be considered as a simulation of noise in the 
inversion procedure having as an input SCIATRAN-generated 
synthetic spectra. 

Results of the inversion procedure described above are 
shown in Fig. 2. It follows that with exclusion of clouds having 
large top altitude and small cloud bottom height (very thick 
clouds), which are rare cases in terrestrial atmosphere, there 
is one-to-one correspondence between retrieved H,. and exact 
He cloud bottom heights. Biases AH = H, — He are given in 
Fig. 3. It follows that biases only weakly influenced by values 
of 7 and generally they are below 0.5 km for clouds having 
the geometrical thickness below 4.5 km. Most of clouds in the 
terrestrial atmosphere have geometrical thicknesses below 1 km 
[3]. Then biases are just 0.25 km as shown in Fig. 3. Therefore, 
we conclude that the technique presented here can be used for 
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Fig. 2. Dependence of the retrieved cloud bottom height on the exact cloud bottom height at r = 10(a),7 = 20(b), 7 = 50(c) for cloud top height 9 km 


(stars), 7 km (triangles), 5 km (circles), 3 km (squares). 


absolute error, km 
absolute error, km 


absolute error, km 
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Fig. 3. 


cloud bottom height (exact), km 


(b) 
Dependence of the absolute error of the retrieved cloud bottom height on the exact cloud bottom height at 7 = 10(a),7 = 20(b), 7 = 50(c) for cloud 


cloud bottom height (exact), km 


(c) 


top height 9 km (stars), 7 km (triangles), 5 km (circles), 3 km (squares). Data are obtained from Fig. 2. 


an accurate estimation of a cloud bottom height from a satellite. 
This is contrary to the general wisdom that satellite techniques 
measure mostly cloud top height properties. 

It follows from Fig. 3 that the accuracy of the retrieval de- 
creases for smaller values of H at the cloud top height equal 
to 9 km. This is related to the fact that analytical equations have 
larger errors for smaller H (see Fig. 1). Note that the probability 
of single clouds having small values of H and large values of 
the cloud top height is low. So values with large biases of the 
retrieved cloud bottom altitude as shown in Fig. 3 do not pro- 
duce significant biases as far as the operational cloud retrievals 
are concerned. 


Itis of importance that the dependence of A H ont is weak. We 
believe that our results are valid for allr > 5, where errors of 
our analytical formulation are small. For smaller thicknesses, the 
approach given here is notneeded because then lidar can penetrate 
a cloud layer and give information on the cloud bottom altitude. 

The results given in Figs. 2 and 3 are obtained allowing for 
the cloud top height measurements uncertainty Ah = +10 m. 
They are not changed appreciably even increasing this uncer- 
tainty in ten times (e.g., to 100 m). The changes of biases as 
shown in Fig. 3 are generally below 30 m then. It means that the 
technique is not very dependent on the vertical resolution of the 
lidar system, which is usually better than 20 m. 
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IHI. CONCLUSION 


We have presented here a new pathway for the determination 
of a cloud bottom height from aircraft or a satellite using simul- 
taneous lidar measurements of the cloud top altitude and the 
cloud reflectance spectrum in the oxygen A-band. Therefore, it 
is of interest to have such a system on a satellite in the future to 
address the climatology of cloud bottom altitudes. Such infor- 
mation on a global scale is needed for climate studies. However, 
it is virtually absent at the moment. 

Most of terrestrial cloudiness is characterized by several 
cloud layers. Our additional numerical experiments (not shown 
here) suggest that the retrieved cloud bottom height is close to 
that of a lowest layer of a multilayered cloud system. 
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